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ABSTRACT 

Spurred by recent observations of 24/im emission within wind-blown bubbles, we study the role that dust can 
play in such environments, and build an approximate model of a particular wind-blown bubble, 'N49.' First, we 
model the observations with a dusty wind-blown bubble, and then ask whether dust could survive within N49 
to its present age (estimated to be 5 x 10^ to 10^ years). We find that dust sputtering and especially dust-gas 
friction would imply relatively short timescales (f ^ 10^ years) for dust survival in the wind-shocked region of 
the bubble. To explain the 24/im emission, we postulate that the grains are replenished within the wind-blown 
bubble by destruction of embedded, dense cloudlets of ISM gas that have been over-run by the expanding wind- 
blown bubble. We calculate the ablation timescales for cloudlets within N49 and find approximate parameters 
for the embedded cloudlets that can replenish the dust; the parameters for the cloudlets are roughly similar to 
those observed in other nebula. Such dust will have an important effect on the bubble: including simple dust 
cooling in a wind-blown bubble model for N49, we find that the luminosity is higher by approximately a factor 
of six at a bubble age of about lO'* years. At ages of 10^ years, the energy contained in the bubble is lower 
by about a factor of eight if dust is included; if dust must be replenished within the bubble, the associated 
accompanying gas mass will also be very important to wind-blown bubble cooling and evolution. While more 
detailed models are certainly called for, this work illustrates the possible strong importance of dust in wind- 
blown bubbles, and is a first step toward models of dusty, wind-blown bubbles. 
Subject headings: ISM:bubbles, dust, stars individual (N49) 



1. INTRODUCTION 

In addition to the energy, mass, and momentum input pro- 
vided by the deaths of massive stars in supernovae, early-type 
stars also launch hypersonic (v ^ 1000- 3000 km s"', whereas 
Cs.iSM ^ 10km s"') winds with mass outflow rates of order 
10"^ to 10"^ yr"' (although clumpin g within the w ind may 



ting the stage for the star's subsequent supernova explosion or 
perhaps gamma-ray burst. 

These wind-blown bubbles and their impact on the in- 
terstel lar medium have bee n observed for many years 
(e.g., I Joh nson & HoggI [19 61 ISmithl \T96% iChu & Laskei 



imply smaller mass-outflow rates; see iHillie^ l2009h . Such 
winds drive forward shocks both into the surrounding inter- 
stellar medium and reverse shocks back into the stellar wind; 
these reverse shocks then heat the stellar wind into a pres- 
surized bubble surrounding the young star, resulting in ob- 
jects known as wind-blown ljubbl e s (W'BBs; for reviews see , 
e.g. iZhekov & Mvasnikovl I2OOOI: IChul I200I lArthuii l2007t 
lChull2008h . WBBs are therefore thought to have four major 
components: (1) an inner, free-flowing, hypersonic wind sur- 
rounded by (2) the hot, post-shocked wind plasma, which is 
enclosed by (3) a thin shell of 10"* K gas which resides in (4) 
the am bient interstellar medium (see Figure 1 of [Weaver et alj 
119771 for a sketch of these components; they also form the ba- 
sis of our Fig.[T4l). 

The structure of these WBBs is useful to study in that they 
help us understand the transformation of gas phases in the 
interstellar medium, and more generally, the impact of mas- 
sive stars on their host molecular clouds, the heating of ISM 
gas, and the input of energy into turbulence in the ISM (e.g., 
[Hensler 2008); they may eve n trigger star formation on their 
periphery (Chu rch well et alJ l2006). Observations of WBBs 
also place constraints on stellar-wind power over the evolu- 
tion of early-type stars. Finally, WBBs are important for set- 
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1980; Naze eLaLl|2001t IChurchwell et al.l l2006l; IWatson etal 
12008). Theoretical work in understanding these ob 
servations has been pursued in paral le l by ma ny re 
(e.g., Johnson & Hogg 1965; 



searchers (e.g., Johnson & Hogg 1965; Mathews 
Piker ner 1968; Dyson 1973; Castor etal. 1975; Falle 
Weaver etal. 1977; Hartquist et al. 1986; Arthur etal. 
. Garcia -Segura & Mac Low 1995; Garcia-Segura et al. 
Arthur et al. 1996; Pittard et al. 2001a b; Naze et al. 



blown bubbles have included detailed modelling of the effect 
of a surrounding inhomogeneous medium (and especially, the 
role of clumps in the surrounding medi um, which will be 
become important later in this work; see [Arthur et"anil993l : 
iPittardet aLllIOOTbl; iNaze etalJl2002l for example). To the 
present day, however, such mod els have not included the role 
of dust. Rec ent observations by IChurchwell et aP (12006 ) and 
IWatson eTa l. (2008) present evidence that dust exists within 
WBBs. What role does this dust play in the structure and the 
evolution of the bubbles, and how might dust modify the im- 
pact of WBBs on the surrounding interstellar medium? 

This paper aims to study the role of dust in wind-blown bub- 
bles in the context of observations of one WBB: N49. We will 
start by briefly outlining the observations of N49 in Section|2l 
This paper will then examine the role of dust within N49; 
in particular, we ask if dust can explain the observed 24/im 
emission around N49 (examined in Section [3]), and then ask 
if dust can survive within the 10^-10^ K, post-shocked wind- 
blown bubble in N49 in SectionH) Then, in Section|5] we offer 
a simple model to explain the observations and our results on 
dust survival. Our conclusions are presented in Section|6l 
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2. OBSERVED MID-IR CONTINUUM DISTRIBUTION 

This paper considers the mid-infrared signature of one 
particular WBB: 'N49.' N49 was disc overed in the in- 
frared and cataloged by Churchwell et alj (12006,) as part of 
a manual search of the Galactic Lega cy Infrared Mid-Plane 
Survey Extraordinaire (GLIMPSE; see lBeniamin et alll2003l: 

turchwell et al. 2009) taken by the Sp itzer Space Telescope 
erner et a l. 2004). Churchwell et all (12006) found over 300 
partial and closed rings in their search of the GLIMPSE data; 
their paper presents an analysis of the characteristics of the 
(incomplete) large population of bubbles reported. A later 
catalog for bubbles within \l\ < 10° in the Galaxy found a 
comparable number. These bubbles were identified by rings 
of 8/im emission that were found to enclose 24 /im emis- 
sion; a few such infrared bubbles had been detected before by 
the Infrared Spa ce Observatory and by the Midcourse Space 
Exveriment (e.g..|Pasauali et al ll2002l ; lDeharveng et al.ll2005l: 
IVasquez et al.ll2009h . 

N49 was l ater th e focus of a more detailed study by 
IWatson et aH (120081) . For that study, N49 was chosen be- 
cause of its relatively symmetric appearance: for instance, 
the 8 /im emission wa s very well fit by a ssuming a spheri- 
cal shell (see Fig. 12 of lWatson et al.ll2008h . This implies that 
perhaps N49 is a useful test case for compa ring with mod- 
els of WBB evolution. 'Watson et al.' (2008) also examined 
24 fim emission observed with the Multiband Imaging Pho- 
tometer aboard Spitzer (MIPS; iRieke et al.l2004 as part of the 
the MIPS Galactic Plane Survey (MIPSGAL; see lCarey etaP 
l2009l) . comparing that data to 8 iim Spitzer obse rvations and 
to 20 cm radio observations (iHelfand et al.ll2006h . These data 
showed that emission at 24 /itm and 20 cm are coincident, shar- 
ing the same annular geometry (although the 20 cm data, in- 
dicating free-free emission from ionized gas, does not drop 
off quite as quickly with radius as the 24 /im emission); both 
were also enclosed by the larger-radius 8 /im emission of the 
photodissociation region (PDR) shell, strongly indicating that 
the dusty post-shocked gas resides within the outlying 8 /im 
PDR shell. However, it is not at all clear if a reasonable mass 
of dust could explain the observed emission, and if dust could 
survive in that post-shocked gas environment. And, if dust 
can exist within the bubble, what role might such dust play in 
the evolution of wind-blown bubbles? These are the questions 
that the present paper sets out to answer 

3. MODELLING THE MID-IR CONTINUUM DISTRIBUTION 

Our first task was to see whether dust emission (alone) is 
a reasonable fit to the observations. To synthesize observa- 
tions of dust in WBBs, we have run photoionization simu- 
lations using 'Clou dy' (version 07.02.01, last described by 
iFerland et al.lll998h . In those Cloudy simulations, we investi- 
gated emission from a constant-temperature ISM plasma with 
an included "normal" dust grain population (using the dust- 
grain size distribution of Mathis et al. 1977) with a dust-mass- 
to-gas-mass rati o of 6.4 x 10"^. We employed Cloudy_3D 
dMorissetl 120061) to map the ID Cloudy simulations into a 
spherically-symmetric 3D bubble that was then projected onto 
the plane of the sky to compare with observations. Note that 
Cloudy cannot model the dynamic post-shock structure of the 
WBB, so for photoionization studies, we model the impact 
of dusty gas on a static bubble of uniform temperature which 
we set. At this early stage, we have no dynamical model of 
dust fluxes into, or out of, the bubble; we are examining the 
emission from a dusty ISM plasma within the bubble, with a 



standard ISM dust-to-gas ratio. 

The parameters for the simulations and analysis are sum- 
marized in Table [T] along with references for the values used. 
We describe and justify those parameter choices in the para- 
graphs that follow. 



Parameter 


Value 


Reference 


central star 


05V 


Watson et al. (2008) 


^wind 


1.5 X lO-^Mgyr-' 


_yink et al. (2001) 


distance 


5.7 kpc 


Churchwell et al. (2006) 


age 


~ 10* yrs 


_Watson et al. (2008) 


^ '^iniier, dusty gas ^ 


19.1 cm-3 


(Best-Fit, 24 ^im data) 


gas & dust density 






power law: r" 


a = -2 


Assumed 


Hnner 


0.6 pc 


(Best-Fit, 24 /xm data) 


''outer 


2.2 pc 


interface of PDR shell 


dust 


ISM grain 


Mathis et al. (1977) 




distiibution 


van Hoof et al. (2004) 


^.bubble 


3.5 X 10* K 


Our dynamical model 






(Sectionim 



Table 1 

The values used for parameters in the Cloudy simulations of the dusty 
wind-blown bubble, N49. 



For the central star and its central st ellar contin- 
uum, we include a star of class 05V (Watso n et alj 
|2008|). Using Cloudy's 'wmbasic' stellar continuum mod- 
els JPauldrach et aTl 12001, which include line-blanketing by 
stellar winds), and calling on the models in Table 1 of 
iMartins et alj (l2005h . we set the star's luminosity to 10^'^' Lq, 
the surface temperature to Tgff = 41540 K, the surface grav- 
ity to \ogig) = 3.92, and the star's metallicity to solar. These 
parameters correspond to a 37.3 Mq star, with a flux of 
hydrogen-ionizing photons of log = lO^'* '* photons cm"^ s"' 
or log Qq = Kf^-^ photons s~' in t o A-k s teradians, very close to 
the value cited bv lWatson et all (12008'). 

The stellar mass loss rate, Mwind is not used directly in the 
Cloudy simulations, but used in our analysis of the velocity of 
the post-shocked gas to estimate (later) the dust-gas friction 
that drags dust out of the post-shocked wind region. For the 
stellar parameter s we'v e adopted, the IDL code 'Mdot.pro' 
from Vink et al. (12001b yields a mass outflow rate of 1 .4 x 
10"^ M0 yr"' and v^o ^ 2600 km s"'. This yields a wind 
kinetic luminosity of 3 x lO"'* ergs s"'. 

The kinematic distance to N49 is 5.7 ± 0.6 kpc 
dChurchwell et al.l 12006). Given this distance, the size of 
the bubble, router, is set by the inner radius of PAH em is- 
sion, so is observationally estimated (i Watson et al.ll2008l) as 
2.2 pc. The density of ambient gas surrounding N49 is not 
well constrained, but unpublished observations of NH3 emis- 
sion (Cyganowski et al., in preparation) surrounding N49 in- 
dicate a high ambient gas density of order 10"* cm"-'. The age 
of N49 is not well constrained, but if the external medium has 
a gas density of order 10"* cm""*, the age of a wind-blown bub- 
ble of the inferred size for N 49 is approximately 5 x 10^ to 
lO*" years ( Watson et aO2008l) . 

Again, since Cloudy does not simulate the heating of the 
wind's post-shocked gas, we start by modelling the WBB with 
a uniform temperature for the post-shocked gas. For most 
of the simulations, we have chosen 3.5 x 10^ K as that con- 
stant temperature; this value for the temperature is informed 
by our own simple dynamical models of wind-blown bubble 
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dynamics with dust cooling (see Section l4~2] i. Changes in the 
WBB temperature affect the amount of dust that is inferred; 
increasing the temperature t o 10^ K as suggested in the dust- 
free models of Freyer et al] (120031) would reduce the amount 
of dust required, but not eliminate it. 

Finally, we assume that the grain distribution, which 
we specify in Cloudy, is an ISM-typ e grain distribution 
jMathis et al.lll977l; Ivan Hoof et al.ll2004l) . We explain below, 
in Section 13.11 why that grain distribution seems to be pre- 
ferred by the observations. 

3.1. Fitting 24 fim observations with Cloudy models 

First, we ask whether dust emission can reproduce the 
Spitzer 24 pm observations of N49. Using the above-defined 
Cloudy and Cloudy _3D models, we simulate the emission in 
the Spitzer MIPS 24 pm bancflby defining a filter centered on 
23.7 /im that spans from 21.65 /im to 26.35 /im; we also de- 
fine, for later use, the MIPS 70 pm band, centered on 7 1 .0 /im, 
spanning wavelengths from 61.5 pm to 80.5 fim. To compare 
to th e data, we take azimuthal av erages of both the Spitzer 
data (IChurchwell et al. 2006; Watson et al.l2008) and the pho- 
toionization model. For the Spitzer MIPS data, we take the 
flux average over 50 concentric annuli that have a Ar of 2 
pixels (with one pixel = 1.25"), centered on the centroid of 
flux within the bubble; the result of this azimuthal averaging 
is shown by the solid line in Figure [T] The background emis- 
sion is estimated by calculating the average flux in an annulus 
with an inner radius of 130 pixels from the bubble centroid, 
and an outer radius of 150 pixels from the bubble centroid, 
where there are no background stars. The la variance in the 
per-pixel flux is shown by the error bars in the figure; plotting 
this variance helps to confirm that the rise (at about 25") and 
fall of the 24 pm emission is present, but with some variations 
that we have averaged out. 

We first performed a large-scale, manual search of param- 
eter space for the approximate gas density and the approx- 
imate WBB inner radius that reproduced the observations; 
we started our search at « = 1 cm"-' and rinner = 0.5 pc (the 
last value is taken from the approximate observed inner ra- 
dius of 24 /im dust emission). After manually increasing the 
gas density, n, to reproduce approximately the observed dust 
emission, we then set up an automated grid search of the pa- 
rameter space stretching from log(rinner/cm) = 18.15 to 18.30 

in steps of 0.05, and with < «dusty gas(nnner) > = 19.00 to 

21.00 cm"^ in steps of 0.25, followed by higher-resolution 
searches with step-sizes down to 0.01 in log(rinner) and 0.0625 
in < «dusty gas('"inner) >■ This allowcd much more accurate 
estimates of the WBB parameters. In all of our parameter 
searches, there was no apparent degeneracy in for the two 
parameters investigated, and only one minimum in was 
found, producing a unique fit to the data (we used Poisson er- 
rors in calculating x^, weighting the large-scale radial decline 
in 24 pm emission with radius, and not the central depression 
in emission). 

We find a reasonable match to the 24 pm emission, as 
shown by the dashed line in Figure [T] For this fit, we spec- 
ify that the dusty gas density scales as n{r) oc r~^ and find the 
best-fit initial density and inner radius (again, see Table [T]for 
input parameters and resultant best-fit parameters). This par- 
ticular power-law for the density was chosen because it fit the 
observed pattern of 24 /im emission with radius; we will offer 

^ http://ssc.spitzer.caltech.edu/mips/descrip.html 



a theoretical explanation of the dropoff in density with a more 
detailed model in Section |5] 

Note that while we include a central hole in the dust dis- 
tribution, the model does not reproduce the strength of 24 /tm 
emission towards the center of N49. This may perhaps be the 
result of some non-sphericity in the N49 bubble, resulting in 
less dust along the line of sight to the central source. Or, as 
we suggest later for other reasons, if the dusty medium inside 
the WBB is clumpy, that could result in an enhanced contrast 
between the central hole and shell of 24 /tm emission for some 
lines of sight (K. Wood 2009, private communication). 
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24 /im Cloudy Model 
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Figure 1. The best-fit model for 24 /im emission within the wind-blown bub- 
ble, using an ISM dust-grain distribution with the photoionization model pa- 
rameters given in Table [T] The solid line shows the azimuthally averaged 
flux from Spitzer MIPS data, while the dashed line shows the photoioniza- 
tion results, using the same azimuthal averaging. The error bars show the Icr 
variation in the flux around the averaged aperture. 

This fit shows that it is possible, with an ISM-like dust dis- 
tribution and with a plausible gas density & temperature, to 
reproduce the 24pm observations. But is an ISM-like dust 
distribution the only one that is consistent? We can check fur- 
ther whether our dust distribution is reasonable by compar- 
ing to observations at other wavelengths; in particular, Spitzer 
MIPS data at 70 jim taken simultaneously with the 24 fim ob- 
servations. If small dust grains are preferentially destroyed, 
then N49 would have a dust-grain distribution weighted to- 
wards large-dust grains; requiring such large grains to repro- 
duce the 24 /im emission would require a large increase in the 
total number of large grains; that increase would produce a 
greater amount of 70 /im emission relative to 24 pm emission 
than for a normal ISM-like dust distribution. 

It is important to point out, however, that the observed 
70/im emission has many contributions: it is a function of 
the total starlight intensity in the vicinity of N49, including 
stars other than the central OV star, and is not due to the 
hot post-shock gas in the WBB, as for the 24pm emission. 
It is therefore not surprising that (as we will see) there is a 
significant background of 70 pm emission observed far from 
the post-shocked region of N49. In addition, the currently 
available 70 /tm Spitzer images have obvious stripping arti- 
facts (i.e., not a background matching problem), one of which 
runs directly through the image of the N49 bubble. However, 
to our knowledge, the 70 pm data have been carefully cali- 
brated and should provide a reasonably correct approximate 
upper limit to constrain the dominance of large dust grains in 
N49. 



4 



Everett & Churchwell 



So, to examine this constraint, we ask whether a grain pop- 
ulation with only large grains could also fit the data. We test 
this idea by finding the best-fit model to the Spitzer lAfim 
MIPS data with a large-grain-only population (defined as only 
having grains of size 0.1 to 1 /^m). For this fit, the predicted 
24 /im emission is very similar to that shown in Figure [T] 
the best-fit initial radius for this model is log(rinner/cm) = 
18.3 and < Wdusty gasC^nner) >= 17.4cm"-'. Most importantly, 
though, to fit the 24 /im data, the dust/gas ratio is lO'* ** times 
greater than the normal ISM dust-grain case, so we are clearly 
forced to an extreme case of a large-dust-grain dominated bub- 
ble to fit the 24 /Ltm emission in a similar fashion as shown in 
Figure [U 

Figure |2] compares the observed 70 /im emission with the 
predictions from Cloudy models with a normal ISM-like dust- 
grain distribution (shown with the dashed line, as in Fig. [T]! 
and the large-grains-only distribution (shown with the dot- 
dashed line). The WBB with large dust grains only clearly 
overpredicts the observed emission (by more than a factor of 
two at the peak of the predicted emission), whereas a WBB 
model with an ISM-like dust-grain distribution is, at its peak, 
a factor of four below the observed emission. In addition, the 
observed 70 emission peaks much further from N49's cen- 
tral star than the 24 /im emission, reaching a maximum very 
near young stars that appear to be situated on N49's rim. 

Returning to the comparison of the fits, this comparison 
shows that, along with the extreme over-abundance of dust re- 
quired to fit the 24 /im emission with a large-grain-only dust 
distribution, the observations seem to require a significant 
population of small dust grains as well. Again, we caution that 
the 70 /im emission contains contributions from all surround- 
ing stellar sources, so we do not expect the Cloudy model to 
match the 70 /im emission, but it is clearly a valuable con- 
straint on the dust population in N49, and other wind-blown 
bubbles (Watson et al., in preparation). 



70 nm MIPS Doto 
70 /im Cloudy Model w/ ISM Grain Distribution 
70 /im Cloudy Model w/ Lorge Groins Only 



- Input Stellar Flux 
Emergent Stellar + Bubble Flux 




40 60 80 100 

Aperture Rddius [Arcsec] 



Figure 2. As in Figure |2] except that data and predictions are at 70 /im, and 
we include a prediction for tlie case wliere dust grains are limited to liave sizes 
of 0. 1 to 1 /im only. Such a large-grain population results in very high 70 /im 
flux compared to the observed flux. The error bars show the Icr variation in 
the flux around the averaged aperture. 



Let us now examine where the emergent 24 /xm flux is com- 
ing from: is it indeed dominated by dust, and how does the 
dust flux compare with the stellar flux and other contribu- 
tions? To examine this, we show the model continua in Fig- 
ure[3] The stellar flux is shown by the solid line (which shows 




0.01 0.1 1 10 
Wavelength [/Am] 

Figure 3. A comparison of the continua from the best-fit photoioniza- 
tion model. The solid line shows the stellar continuum that is input into 
Cloudy. The green, dashed line shows the full continuum emerging from 
the star+wind-blown bubble system. In the X-ray and UV (from 10^ to 
10"- /im), the 3.5 X 10* K gas of the wind-blown bubble dominates, while at 
A > 4 /im, dust emission clearly dominates. 



the input stellar flux); it is reasonably approximated with a 
blackbody, although stellar-wind absorption features (around 
r^ 0.1 /im) are present. The green, dashed line shows the flux 
that emerges from the dusty, wind-blown bubble and star, to- 
gether At wavelengths from 10"^ to 10"^ /im (in the UV to 
X-ray regime of approximately 1 to 100 A), the flux is dom- 
inated by the hot wind-shocked gas in the bubble, where we 
have set r = 3.5 X 10^ K. On the other side of the stellar con- 
tinuum, at long wavelengths of 4 /im and beyond, the dust 
emission dominates the continuum. In fact, one can see the 
10 /im silicate feature superposed on the dust continuum. We 
have checked that no other line features contribute in this re- 
gion, in our model; of the atomic lines in Cloudy, H I recom- 
bination lines in the Spitzer lAfim bandpass contribute more 
than any other atomic lines, and they remain at a luminosity 
nine orders of magnitude smaller than the dust emission. In 
this dusty WBB, dust emission dominates in the 24 /im Spitzer 
bandpass. 

For later comparison with dynamical models, an estimate 
is needed for the amount of dusty gas required in the wind- 
blown bubble to reproduce the observations. Given the pa- 
rameters for the best-fit model, the total mass of dusty gas 
would be 4.8 x lO^-' g, or approximately 2.4 M©. How does 
that compare with the mass available in the area swept up by 
the bubble? For the same outer radius, assuming an ambient 
density of 10"* cm"-', approximately lO'* solar masses would 
be available. Thus, the dusty gas required is a small fraction 
of the gas available. 

4. CAN DUST SURVIVE IN THIS ENVIRONMENT? 

These observations and models present a strong case for 
N49 as a dusty, wind-blown bubble. However, there are a 
variety of processes that could either destroy or remove dust 
from the interior of a wind-blown bubble: sublimation, sput- 
tering, radiative acceleration, and/or gas-dust friction. Can 
dust survive within N49? 

We first examine the possibility of dust sublimation in Fig- 
ure ID In our Cloudy models, the graphite-grain temperatures 
are quite low, spanning the range of 40 to 140 K depend- 
ing on position in the wind-blown bubble and on grain-size 
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(the silicate-grain temperatures are quite similar: the small- 
est silicate grains are approximately 25% cooler at the bub- 
ble innermost radius, and have temperatures within 10% of 
the graphite grains over most of the wind-blown bubble). For 
a given grain size, the variation in grain temperature is only 
about 50% throughout the entire bubble. Overall, the grain 
temperatures are much lower than the grain sublimation tem- 
perature of 1 500 K, so the grains are in little danger of being 
sublimated away. We note that, also, the dust temperatures 
given here are quite similar to those observed in other wind- 
blown bubbles (Watson et al., in preparation). 




Figure 4. The graphite-grain temperature as a function of radius in the wind- 
blown bubble and as a function of dust-grain size. The grain temperatures are 
significantly lower than the grain sublimation temperature of 1500A', so even 
the smallest, hottest grains do not subhmate away. 



We next test whether the grains are sputtered away by gas- 
grain collisions. For this, we first assume that the grains are 
stationary relative to the bulk motion of the post-shock gas 
(i.e., that the gas and dust have the same flow velocities), 
an approximation that we will later see is quite good, given 
the dust-gas coupling in the post-shock region. We calculate 
the spu ttering rates from the formula derived by lTielens et al.l 
(119941) . and present the results in Figure|5] 

The sputtering rates depend only on the size of the grain and 
(for the assumption that the gas and grains are well-coupled 
and moving at the same velocity) the temperature of the sur- 
rounding gas. We can see that for the assumed temperature 
of 3.5 X lO^K, the sputtering times for grains with radii less 
than about 0.1 /im are < 10'' years, the approximate estimated 
age of N49 lO^yrs). We caution that these sputtering 
timescales should be recognized as upper limits on the sput- 
tering time; before dust is accelerated to the gas velocity, the 
dust grains will suffer collisions with gas ions of greater rela- 
tive velocity, and so will have a somewhat lower survival time 
than shown here. 

Of course, this result is strongly dependent on the gas den- 
sity in the post-sho ck region. Integrati ng the equations for 
a simple (dust free) i Weaver et al.l (Il977h wind-blown bubble, 
we find that for an external density of 10"* cm"^, the bubble's 
internal density is approximately lOcm"-^ at an age of ap- 
proximately 5 X 10^ years (we will discuss this model and 
more complicated models later in Section |5]l. This results 
in much fas ter sputtering than in the internal density in the 
I Weaver eTal.. (.1977.) models, where the internal density was 



10=1 

0.0001 



0.0010 0.0100 0.1000 
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Figure 5. The sputtering rates for grains of a variety of sizes within the 
3.5 X 10* K post-shock bubble of N49 for the assumed post-shock gas density 
of 19.1cm"' (the sputtering rate increases linearly with gas density). The 
sputtering times for grains are smaller than the estimated age of the N49 
bubble for grains of size less than 0. 1 ^m. 



approximately 0.03 cm , because of their lower external den- 
sity of 1 cm"-' (at t = 10'' years in their Figure 3). 

We also check how quickly the grains are evacuated from 
the bubble by being swept up in the post-shock gas; the results 
of this calculation are shown in Figure |6] 




1.0 1.5 2.0 

Launching Radius [pc] 

Figure 6. The residence time of dust within the wind-blown bubble, com- 
puted using the frictional force the dust grains feel because of the post-shock 
gas streaming past them. We calculate the time it takes each dust grain to be 
acceleration from its starting position ("launching radius") to the edge of the 
nebula, and plot that time for a range of dust-grain sizes. 



To calculate the acceleration of dust due to dust-gas fric- 
tion, we make the following assumptions. First, we assume 
that the velocit y of the post-shock gas follows the velocity 
law derived bv Weaver et al.l (119771 page 379). As mentioned 
previously, we assume a wind velocity of 2600 km s"'. We 
calculate the friction between dust and gas using the results of 
Berruver& Frisch (1983), and model the acceleration of the 
dust in the presence of that drag force and the gravitational po- 
tential of the central star (we find radiative acceleration only 
leads to a small drift velocity between the dust and gas that 
ranges between 5 and <15kms"', depending on grain size, 
which is negUgible compared to the > 500 km s~' velocities 
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of the post-shock gas and dust, as shown below). We sepa- 
rately derive the acceleration vs. grain size for dust at a range 
of initial radii within the wind-blown bubble, and then calcu- 
late how long it takes a dust grain to transit from its starting 
position to the edge of the nebula. This gives us the grain resi- 
dence time, which is plotted for a range of launching radii and 
dust-grain radii in Figure |6] 

As expected, the larger dust grains have greater inertia, 
and so take longer to be accelerated by dust-gas collisions, 
and have correspondingly longer residence times, but only 
by a factor of approximately two. Again, given the approx- 
imate age of N49 at 5 x 10^ to 10^ years, it appears that 
all dust would be evacuated from the bubble on timescales 
of order 10^ years, much less than the bubble age. In fact, 
gas-grain friction would evacuate dust from the bubble on 
shorter timescales than the sputtering timescale; the sputter- 
ing timescale is smaller only for very small grains with sizes 
of order 0.001 fim. So, both processes indicate a much shorter 
dwell time for dust than the approximate age of the bubble, 
although gas-grain friction is the primary mechanism of dust- 
grain removal for grains of all sizes. 

One possible explanation for the presence of dust in this 
environment (that we will investigate in more detail in Sec- 
tionlU is that dust grains could be the product of evaporating, 
colder, more dense condensations, or 'cloudlets', within the 
wind-blown bubble. Perhaps dusty cloudlets are advectedinto 
the WBB, and are then stripped away by the post-shock gas 
within the bubble. This process of cloudlet advection and de- 
struct ion has, of course, been studied by re searchers before 
(e.g iCowie & McKed [19771 iPittardet all l2001akb: Pittard 
120071) albeit not including the effects of dusty cloudlets, to 
our knowledge. Another possible source of dust replenish- 
ment might be debris disks of low-mass stars formed along 
with the 05V star in N49, which could be ejecting dust parti- 
cles. 

As dust emission seems to fit the observations, we continue 
to pursue the role of dust in cooling the bubble in the follow- 
ing sub-sections. We first show how the charge on dust varies 
within the bubble, and then what effect the dust cooling has on 
the bubble; we will then assemble these results into a model 
ofN49. 

4.1. Charge on Dust in the Wind-Blown Bubble 

The impact of dust on the thermal balance of the dusty 
WBB is determined largely by the charge on the grains. In 
Figures |7] and |8] the average graphite grain charges are pre- 
sented: in Figure|7] the average grain charge is computed over 
all grain sizes and over all radii (with no weighting to take into 
account the greater volume at larger radii, or the greater num- 
ber of smaller-diameter dust particles). Figure |8] shows the 
average charge of grains with radii < 0.01 /im. We investigate 
these trends to better understand the role of dust in WBBs. 

There are two trends to note in Figures |7] and |8] first, as 
the luminosity increases, the increased UV flux ejects more 
photo-electrons from the grains, pushing the grains to more 
positive charges. But, as the average ambient density de- 
creases, the ambient electron density decreases, leading to 
fewer collisions with electrons, and resulting in more posi- 
tively charged grains. Grain charge impacts both heating and 
cooling, which are important in our simulations: collisional 
heating of grains accounts for approximately 32% of the grain 
heating, while radiative heating of grains represents the re- 
maining 68% of grain heating. 

In Figures|2]and[8] we consider only graphite grains, but the 



average charge on silicate grains is very similar to the graphite 
grains; the average silicate-grain charge transitions to negative 
values at slightly higher luminosity and lower density (by ap- 
proximately 0.3 dex), compared to the average graphite grain 
charge, but otherwise the trends are very similar. 
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Figure 7. The average graphite dust-grain charge in our Cloudy models of 
the dusty wind-blown bubble. The average here is taken over all radii, and 
over all grain sizes, and is plotted vs. the luminosity of the central star, L, and 
the density of the post-shocked gas, n/y. (NB: this average is not weighted 
for the volume occupied by grains at each radius, or for the relative popula- 
tions of grains of different sizes!) As expected, at higher luminosities photo- 
emission of electrons becomes an important process and the grains become 
positively charged; on the ordinate, meanwhile, the grain charge becomes 
more negative as fte increases. 
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Figure 8. As in Fig. [7] but with the average taken only over the three smallest 
grain sizes (a < 10"^ /im). The same trends can be seen as in Fig. [7] but for 
the smaller grains, which dominate the grain distribution by number. 



4.2. What is the Effect of Dust on the Wind-Blown Bubble? 

If dust can survive, dust cooling will be dominant in the hot, 
post-shocked gas of a wind-blown bubble. This is illustrated 
in Figure |9] where we compare cooling curves from Cloudy 
where dust has not been included in the model (solid line) and 
where dust has been included (dashed line). We can see from 
Figure |9] that in the hot, post-shocked gas of the WBB (see, 
e.g.. ,Frever et al..2006.) . dust cooling will be dominant for T > 
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Figure 9. Gas cooling curves within tlie WBB for the case where dust is 
not included (solid hne) and when ISM dust is included (dashed line); the 
dust-to-gas mass fraction is set to 6.3 X 10"', as in all of our calculations. As 
expected, where the atomic cooling starts to drop off near 4 X 10^ K, the dust 
cooling starts to dominate. For the hot, post-shocked gas, if dust can survive 
in that region, dust will dominate the coohng. 




Energy [Ryd] 



Figure 11. The ratio of stellar flux transmitted through the wind-blown bub- 
ble to the incident stellai' continuum. The peaks in absorption coiTespond to 
increase in dust absorption properties at 2200 A(~ 0.4 Ryd) and at ~ 1 Ryd. 
The cutoff in the fraction at approximately 20 Ryd is due to the complete 
fall-off of stellar flux at that limit, so the ratio becomes ill-defined. 



increase of mass inside the bubble: 
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Figure 10. The fraction of cooling from dust for a range of values of central 
stellar luminosity and hot-gas density (the parameters for the hot gas in N49, 
inferred from the mass outflow rate, are at the center of this plot). The dom- 
inance of dust is evident across all of the parameter space surveyed; note the 
suppressed zero on the legend, at right. Only on the lower-right portion of the 
plot does Compton cooling start to become a factor 



5 X 10^ K. We can see this for a wider range of parameters in 
Figure [To] 

But how does dust impact the WBB dynamics? First, in 
considering the size of the WBB, the dust column that we find 
is not large enough (in the post-shocked region of the nebula) 
to absorb a significant fraction of the ionizing radiation; over- 
all, dust absorbs less than 3% of the stellar flux in the energy 
range of 0.01 Ryd to 20 Ryd (see Fig.fTTTi. And yet, this small 
amount of dust can significantly impact the cooling and evo- 
lution of the WBB. 

We present a simplified numerical model for N49 by in- 
tegrating the equations of motion for a Weaver et al. (1977) 
WBB and including dust cooling. For this, we use the 
standard eq uations of wind-blown bubble expansion from 
IWeaver eTal . (1977 ) (see, for a summary of the equa- 
tions. i Harper-Clark & Murravll2009l) . with the exception of a 
slightly improved cooling term in the equation for the rate of 



(1) 

where, instead of u sing TBubbie = 2 x 10^ K as assumed in 
IWeaver et al.l (119771) . we use the temperature in the wind- 
blown bubble at each time-step, and calculate the luminosity. 
Lb, using the dust-dominated (at high T) cooling curve from 
Cloudy (note that this equation assumes that the mass-loading 
from evaporation of the surrounding dense shell is dominant; 
we will examine this assumption at the end of Section|5]). To 
start with a simple model, and to isolate the effects of dust 
alone, we assume that dust is added at the same rate as gas to 
the WBB, so that dust is kept at a constant dust-mass-to-gas- 
mass ratio. 

The result of this calculation is showr0 in Figure [12] We 
have checked the c alculation by duplicating the results in 
IWeaver et al.l (Il977h in the dust-free case. Note also the tran- 
sition at f ~ 10^ years, in our models, from dust-free initial 
conditions to the dusty WBB dynamics; the initial radiation 
losses are high enough that Mbubbie ceases to grow in mass, 
and for these times, we set t/Mbubbie/i^f = 0. 

Figure [T3] shows the same model, with the same initial pa- 
rameters, except tha t dust is not include d in this model. One 
can see that (as in IWeaver et alj Il977h . the luminosity in- 
creases with time, in contrast to the model with dust cooling 
presented in Figure[T2] In fact, the dusty WBB has a luminos- 
ity at f = lO** years that is approximately a factor of six higher 
than the luminosity of the dust-free WBB. As a result, by an 
age of 10^ years, if the central star lives that long, the energy 
contained within the dusty WBB is approximately a factor of 
eight less than the energy in the dust-free bubble. 

The temperature indicated by the dynamical model of this 
dusty wind-blown bubble is approximately 3.5 x 10^ K at an 
age of 10^ years; this temperature only slightly decreases over 
time to ~ 3 X 10^ K at 10^ years. This model was therefore 
used to set the fiducial temperature and density that we have 

^ We note that this includes non-radiative cooling of electrons due to the 
electric potential of the dust grains, and as such, we will overestimate the 
observed luminosity of the bubble. 
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Figure 12. Analytic calculation of the growth of a wind-blown bubble with dust cooling included. In the legend, /fBubblc refers to the outermost radius of the 
wind-blown bubble (^Bubble is labeled R2 in Weaver et al.); in N49, this is identified with the innermost radius of the PAH emission. Mflubble shows the amount 
of inass that has been evaporated from the ISM shock and subsuined into the post-shocked gas (Msubblc is Mi, in Weaver et al.). ^Bubble is the total energy in 
hot gas within the post-shocked gas, Rwind is the radius of the wind-dominated, dust-free region inside the post-shocked gas (these are Ei, and Ri in Weaver et 
al., respectively). VBu(,(,le is the velocity of the outermost shock of the bubble, and Lgubblc is the total luminosity (including dust emission and also including 
non-radiative cooling of electrons due to dust-grain charge) from the post-shocked gas (corresponding to V2 and Lj, in Weaver et al.. respectively). The upturns 
in luminosity, energy density, and inass in the bubble at t=200 years are the transitions from Weaver et al. 1 1977) initial conditions to our dynamical model that 
includes dust cooling. The near-constant luminosity is the result of the increasing size of the bubble compensating for both the slowly decreasing density in the 
bubble and the slow drop in cooling as the temperature decreases (this calculation assuines a constant temperature throughout the WBB at any given time, which 
slowly decreases over time). 
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used in earlier sections of the paper, in order to be relatively 
self-consistent with the parameters used. (Earli er models of 
ours assumed a temperature of T ^ lO^K as in Frever et aTl 
(p003); as mentioned earlier, that different temperature does 
not affect our results significantly.) 

It is important to note that, in this model, the gas tem- 
perature inside the dusty WBB is approximately the same as 
the dust-free WBB for ages > lO'' years. This is surprising, 
since dust is such a strong coolant. However, this coolant is 
most important at early times, when the WBB temperature 
is ~ 10^ K; at such high temperatures, the dust cooling rate 
is much higher than the gas-derived cooling rate. So, in the 
early evolution of the WBB, when dust is introduced in the 
model, it is indeed a strong coolant: the resultant WBB lu- 
minosity is very high, and the temperature and energy in the 
bubble drop sharply. In the Weaver et al. model, such a high 
luminosity implies a lower energy input to the outer shell of 
the WBB, and therefore a decrease in the evaporated mass- 
flux from that shell into the bubble. This decrease in mass 
flux and the energy lost due to radiation in this ^ 10^ K phase 
roughly track each other, and so the overall temperature in the 
bubble does not change greatly. Therefore, in this model, the 
dust impacts the growth of the bubble (its size and mass), but 
does not greatly affect the temperature. 

This WBB model also approximately reproduces other fea- 
tures of N49. The inner radius of the bubble reaches 0.5 pc 
at an age of 3 x 10^ years, and at the same time has an outer 
radius of approximately 1.8 pc, close to the 2.2 pc that we de- 
rive from observations. The expansion velocity of the bubble 
is approximately 3 km s"' at that same age. 

But how can dust exist continuously in the post-shocked 
bubble? We consider that possibility in the next section. 

5. A PROPOSED MODEL FOR N49 

Our photoionization and (as yet, quite simplified) dynami- 
cal modeling of N49 shows the possible importance of dust on 
the bubble. However, as shown in earlier sections, due to the 
relatively small timescale of dust advection in the post-shock 
bubble, that dust must be somehow resupplied. As mentioned 
in SectionlH we hypothesize that small, dense ISM 'cloudlets' 
have been over-run by the ionization front of the bubble (see 
Figure \14[ . When these cloudlets interact with the shocked 
wind, they are gradually disrupted, releasing their dust and 
gas into the post-shocked wind. The dust grains from these 
cloudlets then yield the 24 fim emission that is observed. The 
dust adds to the radiative cooling of the post-shock gas, and 
decreases the energy density in the bubble over time (Ep in 
Fig.O relative to dust-free wind-blown bubbles. 

We note that, while the presence and importance of dust 
within wind-blown nebula is a relatively new idea, the pres- 
ence of substructure, or 'clumps', within nebulae has been 
inferred in a variety of other observati ons, from early stud- 
ies of line ratios in the Orion Nebula (lOsterbrock & Flatbed 
Il959h to the impact of clump s on o bservati ons of the radio- 
continuum spectral slopj^ (ICass inelli & HartmaniJ Il977t 
llgnace & Churchwell ^ 2004). Extensive theoretical work has 
already been done on the inclusion of clu mps and the result- 
ing ma s s-load ing (see Pi ttard et al.ll2001a llH): for a review, see 
iPittardI (l2007h and references therein. 

To further check this picture, we estimate how quickly these 
cloudlets are destroyed within the WBB, and check whether 

' The radio-continuum spectral slope can also be modified by the power- 
law drop-off of density with radius tPanagia & FellL1975.;.01non. 1975i) 



such resupply is consistent with observations and the model. 
To examine this, we calculated the mass deposition from the 
cloudlets into the WBB, where the cloudlets are stripped away 
by both ablation and evaporation due to the fast-moving, hot 
post-shocked wind that surrounds them. [The ram pressure 
and thermal pressure of the post-shocked wind would, to- 
gether, suppress photoevaporation of the clouds; in addition, 
even if photoevaporation were present and not suppressed, it 
would result in mass-loss rates of order the ablation and con- 
duction rates, so would not significantly change our estimates, 
especially given other uncertainties (Pittard 2007).] To esti- 
mate the eff ects of these p rocesses, we used the approximate 
equations of lPittardI (12007b . First, for ablation, where the ex- 
ternal flow is subsonic, we calculate: 



ablation ' 



gs 



(2) 



where A4 is the Mach number of the post-shocked wind flow, 
Mc is the mass of the cloudlet, is the speed of sound within 
the cloudlet, p is the mass density of the post-shocked wind, 
and V is the velocity of the post-shocked wind. For a super- 
sonic flow, 



gs 



Meanwhile, the mass loss to conduction is given by 



;2.75 X 10' 



^5/2 - 

7^ gs 



(3) 



(4) 



where oj=l gives the limit of classical evaporation (which we 
adopt here), is the temperature in the wind-blown bubble in 
units of 10^ K, and /"cgc is the radius of the cloudlet in parsecs 
(again, see Pittard 2007|). 

To use these equations, we take the density within the bub- 
ble from our simple dynamical model (shown in Fig. fT2l l at 
f = 5 X 10^ years, and use the approx imation for v within the 
wind-blown bubble from page 379 of iWeaver et all (WfH) : 



^wind 



1 c \ 3/2 p2 
16 



(5) 



again using the wind velocity, post-shocked wind temper- 
ature, and bubble radius from our dynamical model. We 
then check whether there exist reasonable cloud parameters 
(Mc, Cc, and rc pc) that would yield the required mass-outflow 
rate given the above ablation and conduction mass-loss es- 
timates, and which would have enough mass to survive for 
5 X 10^ years. 

Specifically, we require that the cloudlets need to resup- 
ply the observed dusty-gas mass (~ 2.5 M© of gas) every 
10"* years (the shortest timescale for dust sublimation and 
evacuation from the bubble). If the wind-blown bubble is 
'--^ 5 X 10^ yrs old, then 125 Mq of dusty gas is required to 
be resupplied to the bubble interior in that time. The av- 
erage mass loss rate is then, of course, 2.5 x lO'^^MQyr"' 
from all of the cloudlets embedded in the post-shock gas. 
Can cloudlet ablation or conduction (or both together?) in 
the post-shock region match that rate, and with a reasonable 
number of clouds? 

Looking at Equations |2] through H) the conductive and ab- 
lative mass-loss rates depend not only on the number of 
cloudlets, but also the cloudlet density, radius, and temper- 
ature. With these unknown parameters, we cannot claim a 
unique solution. Of all of those parameters, we can only 
approximately set the density within the cloudlets, choosing 
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Figure 13. As in Figure [T2l but without dust cooling included. 
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"cloudlet =10^ cm"^ as an approximate high-density component 
of the external medium. After a search of parameter space, 
we find that requiring these clouds to lose mass at a rate of 
^ 2.5 X 10"^ Mq yr~' can be satisfied if the clouds have a ra- 
dius of approximately 0.05 pc or ^ 10'* AU and a cloud tem- 
perature of about T ^ lOOK. Interestingly, t he size of these 
clouds is i n the range of cloud s observed in iDe Marco et al.l 
( 120061) and lKoenig et all (120081) . The mass of each individual 
cloud is approximately l.SM©; the mass of dusty-gas within 
the bubble could be replenished by a few hundred of these 
kinds of clouds, spread throughout the wind-blown bubble, 
occupying only approximately 0.5% of the volume of N49. 

With these parameters, cloud ablation dominates for most 
of the interior of the wind-blown bubble, whereas thermal 
evaporation due to conduction is important only for R > 
1.5 pc. This transition occurs because the velocity of the 
post-shocked gas decreases with radius (see Eq. |5] from 
IWeaver et al.llT977l) as r ^, so that ablation becomes less im- 
portant with radius in the wind-blown bubble. If ablation 
were the only important cloud destruction process, Mabiation oc 
_y^4/3^,i/3 _ (y/cs ext)'*'^''^'''-' OC v^^^ (where Cs.ext is the sound 
speed in the external medium), and since v oc r"^, we would 
expect Mabiation and the resultant dust density to scale as r~^^^^ 
or approximately as r~^. Conduction at larger radius adds an- 
other mass-loss term, however, so we would expect the dust 
distribution to drop off somewhat more slowly than r~^. Also, 
dust-gas friction will push dust to larger radii within the bub- 
ble and would tend to flatten the radial distribution of dust. 
So, in the end, our model's assumed dust distribution of 
is not unreasonable. Much more work is needed to try to un- 
derstand the details of this process however; we have aimed 
merely to physically justify the assumptions of our model and 
to try to build a reasonably self-consistent picture. It must be 
noted that cloud destruction in this environment is a compli- 
cated and non-linear process, and so the above expressions are 
only rough estimates of this process; the rate of cloud ablation 
would certainly change over time as gas and dust are added to 
the post-shocked gas, for instance. 

One difficulty with the above simple explanation for dust- 
loading that points to further work is the amount of dusty 
gas that is added to the pressuriz ed interior of the bu bble. 
We have employed the models of IWeaver et al.l (Il977h . but 
with 125 M0 of dense cloudlets added to the bubble over 10^ 
years, the mass in the cloudlets will overwhelm the addition 
of mass via thermal conduction from the outer s hell (assumed 
to be the main component of mass addition in IWeaver et al.l 
Il977h . dominate mass addition to the interior of the bubble, 
and cool the bubble interior and impact its dynamics (see,e.g., 
lArdiur et al . 1993; Pittard et al. 2001b; Naze et al. 2002). The 
necessity for so much added gas mass could point to improve- 
ments that must be made to our initial assumptions and calcu- 
lations: perhaps such dense cloudlets have larger dust grains 
than we have assumed that are not destroyed as quickly (but 
that yield smaller dust grains via slow sputtering), or perhaps 
normal-sized ISM grains last longer than the above estimates 
of dust destruction indicate. The possibility of longer-lived 
grains could actually be related to the added gas mass, which 
could cool the bubble interior and impacts the gas dynamics 
within the hot bubble; both of these processes may slow the 
destruction of grains, although, possibly countering this, an 
increased density within the bubble would (by itself) increase 
dust sublimation rates proportionately. 



Regardless, it is clear that the presence of dust and of 
cloudlets point to the necessity of new and more detailed 
models of wind-blown bubbles: when we run simple mod- 
els with the above mass-injection rates (from cloudlets), we 
find that the bubble quickly cools. In fact, even with mass- 
injection rates an order of magnitude lower than those pos- 
tulated above (2.5 x lO~^M0yr"'), the hot-pressure bubble 
cools after only a few hundred years. With mass input rates of 
2.5 X 10"^ M© yr"', the bubble expands for millions of years 
without cooling. This could indicate a critical upper limit for 
the addition of mass in cloudlets. We stress again, however, 
that the above estimates are very simple and point the way to 
necessary improvements of the basic models of wind-blown 
bubble exp ansion, possibly alon g the h nes of models such 
as those of iHanami & Sakashital ( 119871) and iPittardI (120071) . 
or the mode l s of c loudlets within supernova remnants, as in 
iPittard et al.1 dlOOl . 

6. CONCLUSIONS 

We were intrigued by the possibility of dust within the 
wind-blown bubble 'N49' (Watson et al. 2008). Our goal was 
to explain the present observations with dust, and to try to un- 
derstand how dust may survive within that environment. In 
this paper, we found that both the processes of dust-gas fric- 
tion and sputtering of small dust grains can lead to the evacu- 
ation of dust grains from WBBs on relatively short timescales 
of approximately 10"* years. We also found, though, that if 
dust is present, it has a very significant effect on the struc- 
ture of the wind-blown bubble: decreasing the energy within 
the bubble due to dust cooling, and therefore decreasing the 
size of the bubble relative to dust-free bubbles. However, 
in order for dust to be present in N49 today, that dust must 
be replenished; we hypothesize that high-density cloudlets 
(«e 10^ cm"^) can be overrun and subsumed by the wind- 
blown bubble. Those cloudlets are then gradually destroyed 
by ablation and evaporation in the post-shock wind in the bub- 
ble, but on timescales long enough that the cloudlets would 
continue to supply the required density of dust to N49 over at 
least 5 X 10^ years. Such substr ucture has already been sug- 
gested by other researchers (e.g.. lOsterbrock & Flathedl 19591; 
Pittard 2007); what is new here is the consideration of dust. 
We caution that the dust may be supplied by other processes, 
too; for instance, the ejection of dust-grains from the debris 
disks of nearby low-mass Young Stellar Objects might also 
be important. 

This model has several implications. First, the external den- 
sity around N49 must be fairly high to support high-density 
clouds (hc ^ 10^ cm"-'); our dynamical model also assumes 
"exteinai ^ lO'^cm"'' in the external medium. This should be 
tested by constraining the density through ammonia obser- 
vations (Cyganowski et al., in preparation). Also, as dusty 
bubbles age, if the cloudlet density is constant with radius 
within the bubble, dust should gradually be evaporated from 
the inside out, such that there should be wind-blown bub- 
bles with only an outer-rim of 24 /im emission; this could 
be tested by comparing free-free and 24 observations of 
WBBs. Finally, a similar structure of 24fim emission should 
be visible in other wind-blown bubbles; as long the bubble is 
young enough (ages of less than 10^ years for the cloud pa- 
rameters outlined herein) that it has not destroyed all of the 
dusty cloudlets, similar 24 fim emission should be seen; there 
are hints of this already (C. Watson, personal communica- 
tion). Finally, the large amount of mass-loading that appears 



12 



Everett & Churchwell 



embedded""" 
dusty cloudlets 




Figure 14. Our schematic of the dusty, wind-blown bubble of N49. The stellar wind has completely evacuated the central region of dust, leaving a dust-free, 
"wind dominated" zone of radius R„ in the center. We identify the observed central dip in azimuthally-averaged 24 ^m emission in N49 with this wind-dominated 
region. Dust is driven from the post-shock gas by dust-gas friction and destroyed by sputtering on timescales about at least an order of magnitude less than the age 
of the wind-blown bubble; therefore, to resupply dust within the bubble, we hypothesize that small cloudlets are over-run by the expanding bubble, and slowly 
destroyed by heat conduction and ablation. 
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to be necessary to maintain dust in the WBB ma y lead to large 
amounts of cooling and a short lifetime for the IWeaver et alj 
([19771) pressure-driven phase. 
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